Precursor phenomena in frustrated systems 
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To understand the origin of the dynamical transition, be- 
tween high temperature exponential relaxation and low tem- 
perature nonexponential relaxation, that occurs well above 
the static transition in glassy systems, a frustrated spin 
model, with and without disorder, is considered. The model 
has two phase transitions, the lower being a standard spin 
glass transition (in presence of disorder) or fully frustrated 
Ising (in absence of disorder), and the higher being a Potts 
transition. Monte Carlo results clarify that in the model with 
(or without) disorder the precursor phenomena are related to 
the Griffiths (or Potts) transition. The Griffiths transition is 
a vanishing transition which occurs above the Potts transition 
and is present only when disorder is present, while the Potts 
transition which signals the effect due to frustration is always 
present. These results suggest that precursor phenomena in 
frustrated systems are due either to disorder and/or to frus- 
tration, giving a consistent interpretation also for the limiting 
cases of Ising spin glass and of Ising fully frustrated model, 
where also the Potts transition is vanishing. This interpreta- 
tion could play a relevant role in glassy systems beyond the 
spin systems case. 

PACS numbers: 75.50.Nr, 02.70.Lq, 05.70.Fh 



Experiments on glassy systems like structural glasses, 
ionic conductors, supercooled liquids, polymers, colloids 
and spin glasses (SG) Q show that precursor phenom- 
ena occur at some temperature T* well above the static 
transition. In particular, the density-density or spin-spin 
autocorrelation function in glasses M or spin glasses ||] 
has a transition from a high temperature exponential be- 
havior to a low temperature nonexponential behavior at 
T* . Many attempts to relate these dynamical transitions 
to thermodynamic 0] or, alternatively, percolation [^|| 
transitions are present in literature, each one supported 
by numerical simulations [^|-^|. In particular for SG mod- 
els, where random distributed ferromagnetic and antifer- 
romagnetic interactions give rise to frustration, Randcria 
et al. j| suggested that T* should be smaller than the 
Griffiths temperature T c |9[| , i.e. the critical temperature 
of an Ising ferromagnet |1Q] . Intuitively the reason is that 
the randomness of the model allows the presence of ex- 
ponentially rare large unfrustrated regions that reach an 
ordered state at T c , each one with a characteristic length 
and a characteristic relaxation time. The distribution 
of relaxation times below T c gives rise to a nonexponen- 
tial global relaxing correlation for the system, therefore 
T* < T c . Numerically it is found that T* is close to T c 





It has also been suggested that a mechanism which 
would lead to nonexponential relaxation in frustrated 
systems is associated to percolation of the Fortuin- 
Kasteleyn-Coniglio-Klein jl]J (FK-CK) clusters, which 
can be proved occur at a temperature T p < T c . Since it 
is argued that above the percolation transition the avail- 
able phase space become rather ramified, enough to slow 
down the dynamics, this percolation mechanism implies 
that nonexponential relaxation should occur also in fully 
frustrated (FF) Ising systems jL2| where, due to absence 
of disorder, there is no Griffiths phase while the perco- 
lation temperature T p is finite. Recent simulations fTs] ] 
on two dimensional (2D) and 3D FF Ising model in fact 
show the existence of nonexponential relaxation below a 
temperature equal within the numerical precision to T p . 

One might wonder whether also in spin glasses T* = 
T p or instead T* — T c , since both results would satisfy 
Randeira's criterion. Unfortunately it is found that T p is 
less but close to T c and it is very difficult numerically to 
locate precisely T* to distinguish between T p and T c . 
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FIG. 1. Qualitatively phase diagram for the PSG model in 
3D as function of s. Solid lines mark real thermodynamic 
phase transitions, while dotted line marks the vanishing Grif- 
fiths transition. Cross marks Potts vanishing transition in the 
s = I (± J Ising SG) case. An analogous phase diagram holds 
for the PFF model with a FF Ising phase at the place of the 
SG phase [17,18] and without the G phase. 

In order to better understand the role of the percola- 
tion transition and the location of T* in the spin glass 
model in this Letter we consider a frustrated model in 
presence of disorder, where T c is quite larger than T p , so 
that can be easily distinguished. The numerical results 
allow to conclude that T* ~ T c . Moreover, in absence of 
disorder, i.e. in the FF case (without T c ), the relevant 
role of T p in the dynamical transition is confirmed, giv- 
ing a scenario consistent with all the previous results. In 
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particular, the model considered here clarifies the phys- 
ical meaning of T p , relating it to a Potts transition that 
vanishes in the SG case and in the FF Ising case. There- 
fore in the SG case there are two vanishing transitions 
(the Griffiths at T c and the Potts transition at T p ) in 
the paramagnetic phase which give rise to nonexponen- 
tial relaxation. However since T c > T p the first transition 
dominates hiding the effects of the lower transition at T p . 
In the FF Ising model instead there is only the vanishing 
transition at T p which now can be manifested, marking 
the onset of nonexponential relaxation. 

Let us consider the following frustrated model with 2s 
states spin variables whose Hamiltonian is given by: 



H = -sjJ2 [<W 3 - (ajSiSj + 1) - 2] 

(i,o) 



(1) 



where to each lattice site is associated an Ising spin Si — 
±1 and a s-states Potts jl4| spin Oi = 1, ...,s. The 
sum is extended over all nearest neighbor sites, J > is 
the strength of interaction and e,j = ±1 is a quenched 
variable (the sign of the interaction) . 

When the eij are randomly distributed the model is 
a superposition of a ferromagnetic s-states Potts model 
and a ±J Ising SG and we will refer to it as the Potts 
spin glass (PSG); when the configuration is ordered 
in such a way that every elementary cell of the lattice has 
an odd number of -1 signs (i.e. is frustrated) the Potts 
model is superimposed to a FF Ising model, and we will 
refer to it as the Potts FF (PFF) model. For 5 ai<Xj = 1 
(i.e. s = I) we recover, respectively, the ±J Ising SG 
and the FF Ising model. 

Both the PSG and the PFF models exhibit for 
every s > 1 a thermodynamic transition at T p (s) in the 
same universality class of the ferromagnetic s-state Potts 
transition (see Fig.l). This transition corresponds to the 
percolation transition of FK-CK clusters. For example, 
for s — 2 at T p there is a second order phase transition 
in the universality class of the Ising transition 
Therefore T p for s > 1 marks the Potts thermodynamic 
transition, that disappears for s = 1 [jl9| . Both the mod- 
els exhibit also a lower transition in the same universality 
class of the SG transition and of the FF Ising model, re- 
spectively for the PSG |L5) and the PFF model Q. The 
PSG model has also a higher Griffiths temperature at 
T c (s) which is the transition of the model without disor- 
der, namely with all = 1. It is easy to recognize that 
in this case the model corresponds to the ferromagnetic 
2s-Potts model therefore the Griffiths temperature T c (s) 
corresponds to the ferromagnetic 2s-Potts model |Q . 

For s=l, that is the ±J Ising SG, Ogielski (?J showed 
that nonexponential relaxations appear below T c , sepa- 
rating the paramagnetic (PM) phase in a high temper- 
ature PM phase and a low temperature Griffiths (G in 
Fig.l) PM phase. An analogous phase diagram holds for 
the 2D case, with the SG phase suppressed to T = 0. 



Here we will analyze the dynamical behavior of the 2s- 
state PSG and PFF models, with s = 2 in 2D, simulated 
by standard Monte Carlo (MC) spin-flip dynamics. 

PSG model- We have studied the normalized autocor- 
relation function 



/(*) 



Xsc(t) - xsc(t -> oo) 



Xsg(0) - Xso(t -> oo) 
of the time dependent nonlinear susceptibility 



(2) 



Xso(t) 



N 



^2Si(t + to)Si(to) 



(3) 



where N is the total number of spins, io the equilibration 
time, Xsg(0) = N and where the bar stands for the av- 
erage over the disorder and the angular brackets for the 
thermal average. Following Ref. |5) the infinite size be- 
havior of f(t) has been extrapolated at every t plotting 
the data for finite linear sizes Lvsl/L. To test the form 
of f(t) we fitted the data a) with a simple exponential, 
finding good fits only asymptotically for long time and 
for high temperatures, b) with a stretched exponential 
fo • exp[(t/r)' 3 ], finding that it fails to fit the data only 
for short times, and c) with the form f ■ t~ x exp[(t/r) /3 ] 
suggested by Ogielski j^], finding that it fits very well 
the data over all the time's range and the temperature's 
range (see Fig. 2. A) P2"[| . The parameters 8, r and x used 
in the fit b) and c) are plotted in Fig.2.B. In both types 
of fit the correlation functions turn to be nonexponential 
for T < T c . The results shows that T* ~ T c . Since here 
the difference T c — T p is quite large [^3| it is possible to 
exclude T* < T p §|. 

PFF model- We have studied the autocorrelation func- 
tion defined by Eqs. (g, ||) without the average over the 
disorder, since there is no disorder in the model in this 
case. The data, shown in Fig. 3. A, are fitted with the 
above mentioned fitting forms a), b) and c). The fit pa- 
rameters and the integral autocorrelation time pl[ are 
presented in Fig.3.B and show that the onset of non- 
exponential behavior is T* ~ T p within the numerical 
precision J2jj |. 

This dynamical behavior is still present in the s = 1 
case [Q, where the Potts variables are not present any- 
more and the thermodynamic transition at T p disap- 
pears. If there was no frustration the spin variables would 
become critical at T p , like the Potts variables, giving rise 
in the free energy to two minima separated by an infinite 
barrier. However due to frustration in the spin variables 
the two minima in the free energy will evolve into a cor- 
rugated landscape. It is this corrugated landscape which 
gives rise to the nonexponential behavior in the dynam- 
ics. 

In the disordered case one expects that the free energy 
starts to appear corrugated in phase space at the Griffiths 
temperature T c due to the effect of the disorder. As the 
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temperature is lowered at T p the frustration induces more 
roughness in the free energy landscape, however its effect 
will be hard to be detected experimentally or numerically, 
being masked by the effect of the disorder. 

In conclusion, we have compared the dynamical be- 
havior of 2s-states PSG and PFF models, with s = 2. 
These models has two thermodynamic transitions, the 
lower temperature transition being a SG (or FF Ising) 
transition and the higher at T p being an Ising 
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FIG. 2. PSG: A) Nonlinear susceptibility correlation func- 
tion for L — > oo limit from data for square lattices with 
linear sizes (in lattice steps) L = 20, 25, 30, 40 and pe- 
riodic boundary conditions (p.b.c.) at temperature range 
2.25 < T < 3.75 (in J/k B unity), including T c = 3.641 [14] 
and T p = 2.925 ± 0.075 [16]. Since we are well above the SG 
transition (occurring at T — in 2D), relatively low statistics 
produce good data [20]. For clarity we show only some of the 
recorded data for some of the simulated temperatures and an 
enlarged view of the lower temperatures. Points are results of 
the simulation, solid lines the fits with the form c) (see text), 
dashed lines with the form b) and dotted lines with the form 
a). Where not shown, the errors are smaller than the point's 
size. B) Fit parameters used in part A): Circles are the pa- 
rameters for the form c), triangles for the form 6); squares are 
the integral correlation times [21]. Arrows show T p and T c . 



transition. In the disordered (PSG) model nonexponen- 
tial relaxations start to appear at T* ~ T c , i.e. at the the 
Griffiths temperature, excluding the controversial possi- 
bility T* ~ T p thus extending to s=2 the previous 
results for s — 1 |t]]|]- In the PFF model, where there 
is no disorder and the Griffiths phase is not present, the 
effect of frustration which occur cit Tp can be manifested. 
In fact, we find noncxponcntial relaxation starting at T p 
for both s = 2 and s = 1 where a real and, respec- 
tively, a vanishing transition occur. 

This scenario gives a consistent interpretation of all 
the existing results about frustrated spin systems in lit- 
erature, and could be relevant in general to understand 
the high temperature behavior of glassy systems. In 
particular, we suggest that the relation of the precur- 
sor phenomena with the Griffiths phase and the Potts 
vanishing transition should be experimentally studied in 
regular crystals, like cuprates, with doping induced dis- 
order, where both the SG case and the FF case can be 
reproduced. 
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For non integer s the Hamiltonian in Eq. (JlJ) is not de- 
fined, but its partition function, written in a percolation 
approach [A. Coniglio, F. di Liberto, G. Monroy, and 
F. Peruggi, Phys. Rev. B 44, 12605 (1991)], is well de- 
fined for every values of s. In this sense one can consider 
s as a continuous variable and the s = 1 case as the 
limiting case for the s > 1 cases E3[ . 
Data are averaged over a number of random interaction 
configurations n £ [8, 96]. For each temperature and each 
interaction configuration, after 10 5 MC steps (defined as 
the update of the whole system) of equilibration, we have 
recorded the data for a number of MC steps that goes 
from m = 3 • 10 5 for L = 20 to m = 12 ■ 10 4 for L = 40, 
always with n ■ m £ [2 • 10 6 ,6 • 10 6 ], using an annealing 
procedure, i.e. with a slow cooling of the system at each 
temperature. To verify that the statistics is sufficient for 
the considered temperatures we compared the data ob- 
tained starting from a random configuration with those 
obtained from the annealing method, observing the same 
results. 

The integral correlation time is defined as Tint = 
liniw^Ii + E'^ /(*)]• 

In the case c) it is possible to do fits with only two free 
parameters, fitting linearly the short time behavior (de- 
scribed by fo and x) and separately the long time behav- 
ior (described by /3 and r) with fo and x fixed. 
In 2D is T p /(2sa) = l/ln(l + \Msa) with a = 0.0803 ± 
0.003 @ while is T c /{2s) = l/m(l+\/2s) @, therefore 
T c — T p increases with s. 

Same qualitatively results are seen for j3 and x fit param- 
eters associated to autocorrelation functions of energy or 
Potts magnetization, while for these critical quantities r 
exponent diverges near the transition temperature T p . 
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